Interspecific competition underlying mutualistic networks 
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Multiple classes of interactions may exist affecting one another in a given system. For the mutu- 
alistic networks of plants and pollinating animals, it has been known that the degree distribution 
is broad but often deviate from power-law form more significantly for plants than animals. To 
illuminate the origin of such asymmetry, we study a model network in which links are assigned 
under generalized preferential-selection rules between two groups of nodes and find the sensitive 
dependence of the resulting connectivity pattern on the model parameters. The nonlinearity of 
preferential selection can come from interspecific interactions among animals and among plants. 
The model-based analysis of real- world mutualistic networks suggests that a new animal determines 
its partners not only by their abundance but also under the competition with existing animal species, 
which leads to the stretched-exponential degree distributions of plants. 
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Diverse interactions and dependencies among noniden- 
tical elements are characteristic of complex systems Jj- 
Ecological systems are a prototypical example, in 
which numerous species interact via predation, herbivory, 
mutualistic support, competition, cooperation, and so 
on, and their network structure and function have at- 
tracted much attention [3HH]- The mutualistic relation 
between two species is beneficial for the survival and re- 
production of both of them such as animal pollinators 
and flowering plants. Given numerous species of plants 
and insects, producing nectar with different composition 
and flavor and carrying a wide range of preference and 
capacities, respectively, the establishment of individual 
mutualistic relationship depends on specific needs and 
qualification. Nevertheless, the global organization of 
the mutualistic plant-pollinator networks exhibits com- 
mon features [SHTS]. The degrees of plant and animal 
species are distributed broadly in general, but their dis- 
tributions often deviate from a power-law form, more 
significantly for plant species: While the degree distri- 
bution of animals are close to power laws, those of plants 
are of truncated power-law, exponential, or stretched- 
exponential form [TUHl2j . Biological matching, species 
abundance, and the difference between the numbers of 
animals and plants may bring such deviation from scale 
invariance |10( 111) 114) , which, however, remains to be ad- 
dressed further |15j . Interestingly, mutualistic networks 
have their topological features distinguished from trophic 
networks particularly in their nestedness and modularity, 
which is related to the stable architecture varying with 
the type of interaction jl6j . 

To gain insight into the underlying mechanism of mu- 
tualistic community formation, here we propose and 
study a simple growing bipartite network model and ap- 
ply it to analyze real-world mutualistic networks. The 
model is based on a generalized preferential-selection 
rule, being an extension of the model in Ref. [llj . From 
the nestedness and broad degree distributions identified 



in lots of mutualistic communities, the preferential se- 
lection jT7j has been expected to play a role in their 
evolution [TDl [11] [1^1 [Hj . Our model provides a unified 
picture of evolving mutualistic communities. The nonlin- 
ear preferential selection reflects the impact of individual 
characteristics and interspecific interactions on determin- 
ing symbionts. Fitting the model prediction to empirical 
data, we discover a pattern of interspecific interaction af- 
fecting the architecture of mutualistic networks. While 
a new animal species is attracted to a plant species with 
high abundance and thus with many pollinators already, 
it should compete with the existing pollinators. As a 
result, the chance to find a plant species with a large 
number of pollinators is not as high as expected, which 
is shown to cause the degree distribution of plants to take 
a stretched-exponential form. 

To be specific, we consider plant-pollinator mutualistic 
networks as in Fig. [l] (a), in which each node represents 
a species of either animal {A) type or plant (B) type and 
some pairs of nodes of the opposite types are connected 
representing their mutualistic relationship. A growing 
bipartite network (GBN) model defined below and also 
sketched in Fig. [IJb) dH ITMSOj illustrates the topolog- 
ical evolution of the plant-pollinator networks. Initially 
there are £a nodes of type A and £b nodes of type B, 
all pairs of nodes of the opposite types being connected. 
At each time step, a new node of type A (B) arrives 
with probability Pa (Pb), where Pa + Pb = 1- The new 
node of type A [B) are linked to £a {£b) partners of type 
B {A) that are selected with the probability proportional 
to the Xb (Xa) th power of their degrees. Iterating these 
procedures up to time N, one obtains a bipartite network 
of NPa nodes of type A and NPb nodes of type B on 
average. 

The degree-based preferential selection by a new 
species is assumed in this model and can be motivated 
as follows. In case that a new species appears by spe- 
ciation from an existing one, it can be assumed to in- 
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FIG. 1. Plant-pollinator mutualistic network and its model, 
(a) A mutualistic network at a subarctic-alpine site in north- 
ern Sweden [261 130j . (b) GBN model. A new animal (plant) 
node, introduced with probability Pa (Pb) every time step, 
selects Ia {is) plant (animal) nodes one by one, each with 
probability proportional to the As (Xa) th power of its de- 
gree. 



which can be viewed as the effective degree of b for a^""^) . 
On ensemble average, and the selection probability 
will be proportional to k if {r^"^) > 0. On the other hand, 
if {r^'') = 0, (Kb) ^ fc^/''™, in which is the random- 
walk dimension varying with the distribution and cor- 
relation of r''a'^'s '22\. Both cases are considered in our 
model with < \b < 1. Similarly, a new plant b^'^™'> can 
be complementary or substitutable for each pollinator b 
nurturing each animal a, which leads the probability of 
animal a to pollinate 6^"°'^^ to scale with k linearly or 
sublinearly. The impact of such nonlinear preferential 
attachment has been investigated for unipartite growing 
networks [53] and the preferential-selection exponent, A, 
has been measured empirically e.g., for the Internet |24) . 

The mean number Nx(k,T) of nodes of degree k and 
type X evolves with time as [23] 



Nx{k,T + 1) - Nx{k,T) = Px5kj^ + 
Pyiy 
Mx{t) 



[{k - l)^^Nx{k - 1,t) - k^^Nxik, t)] 



(1) 



where Mx{k,T) = E,gv^ fc,^" = Ek=i^ k^"" Nx{k,T). 
We use X and Y to indicate the opposite types of nodes 
such that {X,Y) = {A,B) or {B,A). In the long-time 
limit T 3> 1, Nx{k,T) is proportional to time r such 
that Nx{k,T) — rPxnxik), where nx{k) is the degree 
distribution for type X. It is inserted into Eq. ([!]) to 
yield 



herit the mutualistic interactions of the ancestor species, 
as the protein interactions are inherited by duplicated 
genes 21J. Then, the more partners an existing species 
has, the higher the chance to have a new mutualistic part- 
ner speciated from one of its old partners is. Also, if we 
assume that a new species in a community, appearing by 
either speciation or migration, selects randomly its part- 
ner organisms, the new species will be more likely to form 
a mutualistic relationship with a more abundant species 
due to the existence of more organisms of the species. 
Given that the number of mutualistic partners - degree 
- of a species is positively correlated with the abundance 
of the species, the degree-based preferential selection is 
expected to work in this case, too. 

Restricted resources provided by one type of species 
may draw competition among the other type of species, 
which can give rise to the nonlinear preferential selection. 
For a new animal species a'-"°^' , let us assign ri''-* to each 
pollinator a of each plant &, which can be negative or pos- 
itive depending on the significance of the competitive in- 
teraction between a and a'-'^"^^ relative to the abundance 
of plant b. The sum of those interaction coefRcients over a 
can be a measure of the chance that a^"'^*) can occupy the 
plant b. When the sum is negative, the plant will not be 
selected. Thus, the selection probability of plant b with 
degree k will be proportional to = max{0, X]a=i ^a''"'}: 
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where the constant fix is given by 
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It is worthy to note that Eq. ([T]) has the Kramers-Moyal 
coefficients that are identical at all orders [25] , preventing 
a truncated expansion. While Xa and As are engaged in 
determining /i^ and /i^, respectively, both ha and /i^ 
depend on Pa{= 1 — Pb)t^a, and £b, leading nA[k) and 
nsik) to be coupled. Note that fix = Px/{Py^y) at 
Ax = and Aix = 1 + Pxix/{PYiY) at Ax = 1 (Fig. I 
(a)). 

The dependence of the degree distributions on the pa- 
rameters is manifested in its asymptotic behavior: 
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FIG. 2. Degree distributions in the GBN model, (a) Plot 
of fiA versus Aa with {(.a,(-b) ~ (2,3). (b) Cumulative 
degree distributions Pr(fc_B > k) — "^j, yi^nsik) of plant 
nodes for different values of A_b. {£a,£b) ~ (2,3) for all 
cases, (c) Cumulative degree distributions of animal nodes 
for {Iaj^b) = (2,3) and (1,6) take power laws and k~^^'^ 
asymptotically. A,4 — 1 for both cases, (d) Comparison of 
simulations (points) and theory (lines) from Eqs. ^ and ([3|. 
Here £a = 2,£b = 3, Aa ~ 1, and As = 0.5. For simula- 
tions, we averaged the cumulative degree distributions over 
1000 networks of TV = lO"* nodes. In (a-d), we set Pa = 3/4 
and Pb = 1/4. 
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Remarkably, nx{k) changes from a power-law to a 
stretched-exponential form as Xx deviates from 1 as in 
the example in Fig.[2](b) [23]. For Ax = 1, the degree ex- 
ponent 7x of power-law degree distributions varies con- 
tinuously not only with the fraction of nodes of the two 
types Pa and Pb |11| but also with the initial numbers 
of links £a and £b (Fig-[2jc)). Such sensitive dependence 
in bipartite networks enables us to look into the micro- 
scopic details through the global organization of bipartite 
systems. The numerical solutions are in agreement with 
model simulations as in Fig. [2jd) . 

If £a and £b are integers, the possible value of ■^x 
m Eq. ^ are restricted, which is not the case in real- 
world networks. Considering £x as the mean number of 
initial links, we let a new node of type X assigned [£x\ or 
[£x\ +1 links with probability [£x\ +l-£x or £x - [£x\ , 
respectively, [a; J is the largest integer not larger than x. 
Then Eq. Q remains to be true, with a slight change of 
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kix 



1 



Mx(Li"xJ+l-<?x) 

L<?xJ^x+^x 



,\1X 



n:H.xj,^ (fc>LM), 

{k=Vix\) 
(6) 



with ^Jix = [PxUPy^y)] Y.Zy,,\ k^-nxik). 

We fitted the model predictions of the degree distri- 
butions in Eq. ^ to those of 14 real plant-pollinator 
networks by adjusting £a,£b,Xa, and Xb through 
the combined use of the maximum-likelihood estimation 
and the Kolmogorov-Smirnov statistic ^27j. The fitted 
degree distributions are in good agreement with data 
as shown in Fig. [3] (a). The initial number of links of 
plants £b turns out to be much larger than that of ani- 
mals £a probably due to there being more animal than 
plant species; their normalized values, £aNPa/L and 
£bN Pb / L, are distributed uniformly as shown in Fig. [3] 
(c). Most importantly, the preferential-selection expo- 
nents show a significant asymmetry: Xb < Xa in most 
networks as shown in Fig. |3] (d) . This implies the signifi- 
cantly strong competition between a new animal and ex- 
isting pollinators contrary to the relatively weak competi- 
tion between plants. Thus one can see that the restriction 
on the number of available plant species is a more cru- 
cial factor in shaping the mutualistic community than 
the restriction on available animal species, possibly re- 
lated to the difference in their survival and reproduction 
rates. Consequently, plants with large degrees have the 
advantage of their high abundance screened by the com- 
petition between animals characterized by As less than 1, 
which leads the degree distribution to take a stretched- 
exponential form nB{k) ~ g-(const.)/c b ^ -^^^ remark 
that these observations are not the case for "Hocking" in 
Fig. |3] (d) , for which it has been reported in Ref. [28] that 
the competition between plants is more significant than 
between pollinators, implying Xa < Xb ■ 

To check further the robustness of our results, we ana- 
lyzed the product-country networks representing which 
country exports which product [29]. Both the preva- 
lence of a product and the wealth of each country depend 
on trade volume and thus a country (A) and a product 
(B) may have mutualistic relation in the global economic 
ecosystems. Although Pa < Pb and £a > £b contrasted 
to plant-pollinator networks, we find again that Xb < Xa 
suggesting the same patterns of competition among coun- 
tries as among pollinators. 

In conclusion, we found that the connectivity pattern 
of bipartite networks is sensitively dependent on model 
details. Such sensitivity can be of potential use in con- 
trolling transport and epidemic affected by the connectiv- 
ity in multiconnected systems and, in the present work, 
allowed us to illuminate the interaction patterns hidden 
under the bipartite architecture of mutualistic networks. 
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FIG. 3. Fitting model degree distributions in Eq. (|6| to real-world networks, (a) Degree distributions of plant-pollinator 
networks from different communities (points), denoted by ID made after the first author of each reference [26|, and their 
maximum-likelihood fits (lines) of the GBN model. There are 8f to 1881 nodes in these networks, (b) Degree distributions 
of the product-country networks of 1621 products and 132 countries in different years. If the export volume of a product 
by a country is in the top 25% among all the export volumes, they are connected, (c) Normalized initial numbers of links. 
They are rather uniformly distributed mostly in the region Ia^PaIL + IbNPb/L > 1. (inset) For plant-pollinator networks, 
(■A < £b except for two networks, "Kato" and "McMuUen" . For product-country networks, £a > £b- (d) Preferential-selection 
exponents. Xb < Xa in most networks. For Hocking plant-pollinator network, the competition between flowers is known to be 
more significant than between pollinators [2^, leading to Xa < Xb in our model. 



The stretched-exponential form of the degree distribution 
of plants is shown to be driven by the significant compe- 
tition between pollinating animals, which is true also for 
countries in economic mutualistic networks. Our results 
demonstrate the importance of taking into account the 
cross-talk between different types of links interwoven in a 
system in order to understand and predict its properties. 
While we focused on the global network architecture of 
mutualistic communities, it is highly desirable to extend 
our model to account for the population dynamics and 
diverse interspecific interactions. 
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